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We study propagation of pulses along one-way coupled 
map lattices, which originate from the transition between two 
superstable states of the local map. The velocity of the pulses 
exhibits a staircase-like behaviour as the coupling parameter 
is varied. For a piece- wise linear local map, we prove that the 
velocity of the wave has a Devil's staircase dependence on the 
coupling parameter. A wave travelling with rational velocity 
is found to be stable to parametric perturbations in a manner 
akin to rational mode-locking for circle maps. We provide 
evidence that mode-locking is also present for a broader range 
of maps and couplings. 



I. INTRODUCTION 

An interesting feature of coupled map lattices (CML) 
jjJU is the widespread occurrence of the so-called spatio- 
temporal chaos (^^] , which is irregular behaviour in time 
as well as space. Equally interesting is the appearance 
of coherent structures from an apparently decorrelated 
medium [Q . For example, an interface separating two 
different phases may travel along the lattice ||- |T]]] . The 
movement of such a front depends on the strength of the 
coupling between lattice sites. This paper is concerned 
with investigating some aspects of this phenomenon. 

We consider an infinite collection of cells with a local 
dynamical variable Xt(i), characterizing the state of the 
i-th cell at time t. Both i and t are integer. The global 
state of the lattice at time t is given by the state vector 

X t = {x t (i)} = (..., x t (i - l),x t (i),x t (i + 1), . . .). (1) 

The state of the i-th cell at time t + 1 given by 

x t+ i(i) = Fi(. . .,x t (i - l),x t {i),x t (i + 1), .. .), (2) 

where the function Fi is determined by the local dynam- 
ics in each cell as well as the interaction between cells. 
A physically meaningful interaction will have typically a 



*e-mail: R.Carretero@ucl.ac.uk 

^New address: Center for Nonlinear Dynamics and its Appli- 
cations (CNDA), Dept. Civil and Environmental Engineering, 
University College London, Cower Street, LONDON WC1E 
6BT 



limited range, with decreasing strength for distant neigh- 
bours. Equation (|^) defines the components of the global 
map F: 

F=(...,F_ U F ,F 1 ,...) X t+1 =F(X t ). 

A mapping of global states that will be of use in the 
following is the shift mapping 



G({x t (i)}) = {x t (i + l)}. 



(3) 



It should be clear that F and G commute. 

We consider homogeneous CML (same map Fi for all 
cells) with linear nearest-neighbours interaction. The two 
most widely used models are 

x t+1 (i) = (1 - e)f(x t (i)) + £ - (f(x t (i - 1)) + f(x t (i + 1))) 

(4) 



and 



x t+1 (i) = (1 - e)f{x t (i)) + ef(x t (i - 1)), 



(5) 



which are called diffusive and one-way CML respectively. 
The coupling strength e (also called the diffusive param- 
eter) satisfies < e < 1 to ensure the conservation of 
flux between neighbouring sites. 

In this paper we consider the propagation of localized 
wavefronts along a linear one-way CML. Localized states 
and step states (localized wavefronts) are introduced in 
the next section, where we establish some of their basic 
properties. In section 3 we introduce a one-parameter 
family of piece- wise linear local maps, for which wave- 
fronts can propagate coherently in a simple way. In the 
following section we show that the dynamics of an im- 
portant class of wavefronts can be reduced to the study 
of a discontinuous piecewise- linear circle map, depending 
on two parameters. We study this family with methods 
of symbolic dynamics. These results are then applied to 
the pulse propagation in section 5. We show that for 
each value of the parameters the velocity of the pulse is 
well-defined, and we use it to characterize the structure 
of the parameter space. The latter is found to be struc- 
turally similar to that of a subcritical circle map, with 
velocity playing the role of rotation number. We estab- 
lish the occurrence of a critical transition from a stable 
standing signal to a stable travelling one, as the coupling 
parameter is increased, and we prove parametric stability 
of pulses travelling with rational velocity. The velocity of 
the travelling pulse is found to have Devil's staircase- like 
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dependence on the coupling parameter (see figure 1). In 
the last section we indicate how to generalize our results 
to a broader range of functions and couplings. 



II. LOCALIZED STATES AND PROPAGATING 
FRONTS 

We first define localized states, step states and their 
velocity of propagation and then establish some of their 
basic properties. In what follows, we deal with a linear 
one-way CML, unless stated otherwise. 



A. Localized states and step states 

Let Ax t (i) = \x t (i + 1) - x t (i)\. The mass M(X t ) of 
a global state X t is defined as the variation of the local 
states 



M t = M{X t ) 
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Ax t {i). 



We restrict our attention to states with positive finite 
mass (a zero- mass state is a uniform state). 

For each lattice site i £ Z we define the probability 



Pt(i) 



Ax t (i) 
M t 



(0) 



The mean fi and the variance a of the variable i with 
respect to the probability distribution (^|) are given by 

oo 

fit = K x t) = Yl ipt ^ 

i— — oo 

o\=o{Xtf= £ {i-ii t f Vt {i). 

i— — oo 

The quantity fit is called centre of mass of the state Xt , 
and at its width. The centre of mass and the width give, 
respectively, the average position and spread of the lo- 
calized state. 

A state X t with finite centre of mass and width is said 
to be localized. If, in addition, there exists a positive 
constant k < 1 such that Axt(i) — 0(k'*'), then the state 
is exponentially localized. For an exponentially localized 
state, not only the mean and variance exist, but also all 
central moments of the distribution @. 

We are interested in localized states that can model 
wavefronts. So we define a step state to be a localized 
state for which the configuration {x t (i)} is asymptotic, 
for i — > ±oo, to two distinct fixed points of / fl2|: 



The mass of a step state is bounded from below 

M(X t ) > Ax* = \x* + -x*_\. (8) 

The simplest such state is a pure step state, that is, a 
step function between the two fixed points 



P(k) = x t (i) = 



if i < k 
if i > k. 



(9) 



lim Xt(i) 



/04 



(7) 



For this pure state /x(P(fc)) = k and a(P(k)) = 0. 

A step state has minimal mass or minimal variation, 
i.e. M(xt) — Ax* (cf. (||)), if and only if its local states 
are ordered, that is if Xt(i) < Xt(i + 1) when x*_ < x* + (or 
x t (i) > x t (i+l) when x*_ > x+) for all i. It is plain that if 
the configuration of the one-way CML has minimal mass 
at time t it is going to preserve its minimal mass at any 
time greater than t provided that / is a non-decreasing 
function if x*_ < x* + and non- increasing if x*_ > x* + . 

The global dynamics causes centre of mass and width 
of a localized state to evolve. In this respect, two points 
need consideration. Firstly, the existence of the proba- 
bility pt does not automatically imply that of pt+i, and 
secondly, the image of a localized state is not necessarily 
localized. To achieve the above conditions, we impose 
some mild restrictions on the local map / as well as the 
state X t . The importance of the states defined above is 
that they are invariant under the dynamics of the CML 
— for the complete statements and proofs see appendix 
A. This gives the framework for the existence of steady 
step states and travelling ones. 



B. Travelling velocity and basic properties 

From the previous definitions and results it is clear 
that if the initial state Xq is a step state, then fi(Xt) and 
c(X t ) are defined for all t > (see appendix A for de- 
tails). Our main interest is to determine the average ve- 
locity v of the centre of mass along a global orbit, namely 

v = v(e,X ) = lim \(fi t - Mo). 

t^oo t 

A formal analysis on the stability of exponentially local- 
ized step states establishes the following quite general 
property of the velocity of a step state in one-way CML: 



Theorem 1. Let Xo be an exponentially localized step 
state of a one-way CML. Let the local map f be bounded, 
and let f be a contraction mapping in a neighbourhood 
of the fixed points (0). Then, for all sufficiently small e, 
v(e) = and v(l — e) = 1, independently of Xq. 

The proof of this theorem is given in appendix A. 
The condition that / be a contraction mapping is a bit 
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stronger than the condition that x± be attracting. We 
shall see that for an initial localized interface of the type 
, the critical behaviour is much richer than a discon- 
tinuous transition from v = Q to v = 1, at some inter- 
mediate value of e. Rather, the dependence of v on e is 
characterized by infinitely many critical values, in corre- 
spondence to the boundaries of intervals in e, where the 
velocity is a given rational number (see figure 1). 

A general feature of the velocity of the travelling wave 
is its symmetry. Let us define 5 = 1 — e and consider a 
moving reference frame with a unit positive velocity: 

{y t (i) = x t {i) 
y t +i(i) = x t+ i(i + 1) 
Vt+k{i) = x t +k(i + k). 
The one-way CML in the moving frame now reads 

y t+ i(i-l) = 6f(y t (i)) + (l-6)f(y t (i-l)) 
=>. y t+ i(i) = (l-6)f(y t (i)) + Sf(y t (i + l)). 

(10) 

Equation ([l0|) represents a one-way CML of the type (||) 
but coupled in the opposite direction. Thus, taking into 
account the moving reference frame and the definition of 
S we can conclude that 

1 - v (e) = v(l - e). (11) 

Therefore the velocity of the travelling wave in a one- 
way CML of the form (j|) is symmetric with respect to the 
point (e,v(e)) = (1/2,1/2) (see figure 1). This important 
symmetric property allows one to restrict the study of the 
velocity to the e-interval [0, 1/2]. 



III. A PIECE- WISE LINEAR LOCAL MAP 

In this section we deal with the coherent propagation 
of a wavefront (without damping or dispersion) along the 
lattice. We confine our attention to states with minimal 
mass. In order to ensure that this property is preserved 
under iteration, we shall require the local map / to be 
non-decreasing as stated before. In addition, we ask / to 
be continuous and to have two stable fixed points x*_ and 
x* + (with x*_ < x*l_). So / will have precisely one unstable 
fixed point at x* with x*_ < x* < x* + ■ The local dynamics 
takes place on the interval / = [a;^_,a;i_], where there are 
two basins of attraction I_ = [x*_, x*) and I + = [x* , x*,] 
for the fixed points x*_ and x+, respectively. 



A. The local mapping 

A simple one-parameter family of maps that fulfills the 
above requirement is given by (see figure 2) 

{-1 if x < -a 
\x ii -a<x <a <a< 1, (12) 
1 if x > a 

and by the step function between —1 and 1 centered at 
the origin when a = 0. The fixed points x*_ = — 1 and 
x* + = 1 are superstable, with basins 7_ = [—1,0) and 
/+ = (0, 1], respectively. The repeller is the origin x* — 
0. The coupled map lattice now depends on the two 
parameters e and a, and the parameter space is the unit 
square. The one-way CML with local mapping f a and 
coupling e will be denoted by F ea . 

We partition the interval I = [—1,1] into the unstable 
domain U = (—a, a) and the superstable domains S- — 
[— 1,— a] and = [a, 1]. Any site falling within S± 
maps, under / a , to x± at the next iteration. 

We begin considering states of minimal mass Xt whose 
configurations {xt(i)} form a non- decreasing sequence. 
Because X t is a localized state, only a finite number n of 
local states belong to U : 

ie{k+l,...,k + n} x t (i)eU. (13) 

For i < k (i > k + n), the local states belong to the 
superstable domain S- (S*+) and their image under f a is 
— 1 (+1). Consequently, if we let X' t be the state obtained 
from X t by replacing Xt(i) with —1 for i < k and with 
+1 for i > k + n, then F e . a {X t ) = F e>a {X' t ). The state 
X' t is the reduced state associated with X t , and we write 
X t ~ X' t . We will reduce all states after each iteration of 
the mapping F e ^ a . 

So, without loss of generality, we shall restrict our at- 
tention to states of the type 

X t = {...,-l,-l,x t {k + l),...,x t (k + n),l,l,...) 
with Xt(i) < a, i = k + 1, . . . , k + n. 

The number n of local states in U is called the size of 
X t . The knowledge of the range (|l3| ) provides useful 
information about position and width of a step state. 
Indeed the center of mass of X t satisfies the bounds 

nil — a) , \ nil — a) 

k+ K 2 > <»{X t )<(k + n)- K 2 1 (14) 

and the width 

nil — a) .„ . n . . 

- 2 ^ < <X t ) < -. (15) 

The image of X t is 

x t+ i = (...,-i,-i,/-Mfc + i)),..., 

x t+ i(k + n), f+(x t (k + n)), !,-■■) 



(I) 
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where 



f+(x)= |z + (l- e ). 



(16) 



X t+ i has at most n + 1 sites in £/. The reduced state 
X' t+1 is one of the following: 



(a) (...,-1,-1, x t+1 (k + l), x t+ i(k + 2), 
x t+ i(k + n), 1, !,•••) 



(6) (■••,-1,-1, 



-1, x t+ i(k + 2), 



x t+ i(k + n), x t+ i(k + n+ 1), l,---) 

(c) (...,-1,-1, x t+1 (fc + l), x t+ i(fc + 2), 
x t+ i(fc + n), x t+ i(fc + n + 1), 1, ...) 



(17) 



Loosely speaking, as the time increases from t to t + 1 , 
the state has not propagated in case (a) , it has advanced 
to the right by one site in case (6) and it has propagated 
with dispersion in case (c). 

To decide among these possibilities, one needs to in- 
vestigate the value of /_ and /+ at the boundary sites 
i = k + 1 and i = k + n, respectively. Letting 



7- 



a(e 



1 



7+ = 



a(a + e — 1) 



one verifies that the functions f± satisfies the following 
equations 



f-(-o) = -1 
/_(o) = l-2e 
!-{!-) = -a 



/+(-«) = 1 
/+(«) = 1 
/+(7+) = a. 



2c 



(18) 



We define the gap of the CML with local map f a to be 
the open interval T = (j + , 7- ) . Its length (with sign) is 
given by 

, . a(l - a - 2e(l - e)) 
7 (e, a) = 7 _ - 7+ = ^ U-. 19 

e(l - e) 

Let x eU. Then f+(x) € U if x < 7+, and /_(x) € C/ 
if a; > 7-. From this and the fact that Xt(i) is a non- 
decreasing function of i we obtain that the size of a step 
state of minimal mass cannot increase, unless it is zero, 
for a gap of non-negative size. 

In particular, it is not difficult to prove that when the 
gap size is non-negative, if the size is not greater than 1 
at time 0, then it will remain so for all times. A minimal 
state is a minimal mass step state of size or 1. There- 
fore, if we start with a minimal state, the configuration 
of the CML will remain a minimal state at any time. 



IV. SYMBOLIC DYNAMICS OF MINIMAL 
STATES 

We develop a symbolic description of the dynamics of 
minimal states in the case of non-negative gap length. 
Their dynamics can be reduced to the iteration of a one- 
dimensional piecewise linear circle map — the auxiliary 
map. 



A. The auxiliary map 

Without loss of generality, we consider minimal states 
of the form 



X t = (...,-1,-1,^^,1,1,...), 



(20) 



where either |xt(i)| < a or Xt(i) = —1, for which we 
introduce the shorthand notation 

X t = [x,i] t , 

(the subscript t will often be dropped). Thus [— = 
P(i) — cf. @. The image of X t is given by 

X t+ i = (.-., -1,-1, f-{x t (j)), f+(x t (i)), 1, • ■ •)• (21) 
We define 

1 - a 



(22) 



and note that e c is the area enclosed between the graph 
of f a and the diagonal, for x* < x < x* + . The following 
two results characterize completely the evolution of X t 
in the parameter ranges e < e c and e > 1 — s c . 



Theorem 2. Let X t be a minimal state of the form (|20]), 
and let 



t(i 



1 — e — a 



a(l-e) 



a — e 



(23) 



(i) If e < £ C) then for x(i) < x_ (x(i) > X-) the state 
[x,i] reaches P(i) (P(i — 1)) in a finite time. The state 
[x—, i] is fixed and unstable. 

(ii) If £ > 1 — e c , then for x(i) < x+ (x(i) > x+ ) the 
state [x, i) reaches P(i + k) (P(i + k — 1)) in a finite time. 
The state [x+ , i] is Gxed and unstable under Go F, where 
G is the shift mapping (||). 

The proof of the above results is given in appendix B. 
From Theorem 2 it follows that the velocity of the trav- 
elling interface satisfies 



v(e) 



if < £ < £ c 

1 if l-£ r <£<l. 



(24) 



4 



As a consequence, we have that if e < e c then F(P(i)) ~ 
P(i), and if e > 1 - e c then F(P(i)) ~ P(i + 1). In the 
rest of this chapter we shall assume that e c < e < 1 — e c . 
We partition U into three intervals, namely 

U+ = {x : —a < x < 7 + } 
r = {x : 7+ < x < 7_} 
£/_ = {x : 7_ < x < a}. 

Because we are assuming the gap length to be non- 
negative, we are dealing with a minimal state and then 
the three possibilities for the state ( ^l|) are: 

(a) x t e U- f-(x t ) e U, f+{x t ) e S+ 

(&') i,er f_(xt) G SI, G 5+ (25) 

(6) x t ec/+ f_(xt) e Ufa) e U. 

corresponding to the reduced states — see ( |l7|) 

(a) X t+ i = (..., -1,-1,^+1, 1, 1,...) 
(6') X t+1 = (...,-1,-1, -1, 1, 1,...) 

• • i -1, -1, -1, xt+i, !,•••) 



(6) X t+1 = (. . 
For the case (&') we have the evolution 

[x,i] t — > [-l,i] t +i — ► [1 - 2e,i + l] t+2 . 



(26) 



All possibilities are accounted for by defining the auxi- 
liary map ^ et a to be the 2-parameter map (see figure 3) 

(f+(x) ifx€U + 
$ s ,a( x ) = \ !- 2e iixeT e c <e<l-e c . (27) 

[ /_(») if a; G li- 
lt is plain that $ ea maps ?7 onto itself and if, in addition, 
7 is positive, then U-, T and U+ have positive measure. 
This characterizes the domain and range of <fr e , a - 

The auxiliary map reduces the dimensionality of our 
system. Originally the CML is an infinite dimensional 
system since it has an infinite number of sites. But in 
the case of minimal states the dynamics of the whole 
lattice is reduced to that of the auxiliary map <fr £ . a , as 
illustrated schematically in figure 4. Iterating the CML q 
times amounts to applying the auxiliary map q times on 
the interface site. During these q iterations, whenever the 
branch /_ is used the interface is not shifted; whenever 
f + is applied we shift the interface by one site to the 
right; and finally, when the interface site falls into T the 
CML essentially undergoes two iterations — cf. (26) - 
but the final result is that the interface site value is 1 — 2e 
as well as being shifted one site to the right. 



B. Symbolic dynamics 

The following binary symbolic dynamics for the auxil- 
iary mapping <& £ja will play a decisive role in the rest of 
the paper. 

We assign the symbols '0' to U-, the symbol '1' to U+, 
and the symbol '10' to L. This choice originates from the 
prescriptions (|2^), according to which the location of a 
minimal state remains unchanged for x G J7_, increases 
by one for x S U+, and increases by one for every two it- 
erations for x S L. To the orbit of the auxiliary mapping 
with initial condition x we associate the semi-infinite se- 
quence of symbols 

S = S(x) = (si,s 2 ,s 3 ,...) s fe e {0, 1} 

where we stipulate that the symbol '10' generates two 
binary symbols: '10' — >'1','0'. For this reason, our code 
will not be unique. The space of binary sequences is 
denoted by £2, and is equipped with the usual topology. 

Any orbit with an element in T will be called a gap 
orbit. If this orbit is periodic, then the period of its 
symbolic representation is one greater than that of the 
orbit, because the gap T corresponds to two symbols. 

The rotation number |l^ , p^| of the orbit through x is 
defined as 



1 T 
lim — > s t 

t=0 



(28) 



provided that the limit exists. If this is the case, it is 
plain that < v < 1. 

We will show that the symbol sequence S(x) can be 
characterized in terms the spectrum of i>, which is defined 
as the following sequence of integers 



Spec{v) = \y\ , \ 2v\ , [3^J , . . 



Oil) ^2, 0,3, 



(29) 



where ['J is the floor function. Because < v < 1, 
consecutive elements of Spec(v) differ by zero or one. 
Thus defining 

r(oi, as, . . . ) = (02 - ax, a 3 - a%, ■ ■ ■ ), 

we find that T(Spec(v)) is a binary sequence. 

Using symbolic dynamics it is easy to follow the evolu- 
tion of the minimal state: every iteration corresponds to 
a new symbol in the sequence. If the symbol is '0' then 
the minimal state remains in the same place and if the 
symbol is '1' the minimal state is shifted one site to the 
right. We can assume the velocity of the travelling front 
to be constant, so that the spectrum of the velocity gives 
the position of the minimal state at every iteration, since 
the location of the minimal state can only take integer 
values. Thus r applied on Spec{v) gives a when the 
position of the minimal state remains unchanged and a 1 
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when it jumps from on site to its right neighbour. There- 
fore the symbolic binary sequence associated to a given 
rotation number v can be computed as 

S = t (Spec(v)) . (30) 

If v is rational, then T(Spec(v)) is periodic, as easily 
verified. For example, the sequence S(3/7) is given by 

5(3/7) = t (5pec(3/7)) 

= r(0,0,l,l,2,2,3,3, ...) 
= (0,1,0,1,0,1,0), 

where the bar denotes periodicity. 

The symbolic dynamics of $ £ja is that of a uniform 
rotation with rotation number v analogous to that of a 
subcritical circle map, as we shall see later. The dy- 
namics of the whole lattice is then reduced to the study 
of the one-dimensional map <& ea whose rotation number 
gives the velocity of propagation of the travelling front. 
Using arguments similar to that given in section 1.14 of 
||l5f , it is possible to prove that the rotation number is a 
well-defined, non-decreasing, continuous function of the 
parameters (e, a) and it is independent of the initial con- 
dition. The difference here is that the auxiliary map is 
a non-decreasing function and not a strictly increasing 
one as in Jl5| , but the same line of proof may be fol- 
lowed. Therefore, when the gap size is non-negative, the 
velocity v(e) of the travelling interface is a well defined, 
non-decreasing, continuous function and does not depend 
on the choice of the initial minimal state. 



V. STRUCTURE OF PARAMETER SPACE 

We now turn to the study of the parameter space of 
the coupled map lattice F e ^ a , in the case where the gap 
size 7(e,a) is non-negative (cf. (|l9|)). Then the gap orbit 
is defined, and we denote by v — v(e,a) the velocity of 
the center of mass of the corresponding minimal state. 



A. Mode-Locking 

Using the symbolic coding we now prove that the gap 
r forces the velocity of the travelling wave to be locked 
to rational values if & E ,a has an eventually periodic or- 
bit containing the point 1 — 2e. Mode- locking guarantees 
stability of the propagating signal with respect to para- 
metric perturbations. 

The symbolic coding of the orbit of the auxiliary map 
®E.a gives the velocity of the travelling interface via the 
rotation number (p8[). The presence of a non- negative 
gap r induces a mode- locking of the rotation number in 



the following way. Let us consider the orbit generated 
from an arbitrary point xq in T. If after q iterations x q 
falls into T, the orbit is eventually periodic [jl5| of period 
q with rotation number p/q, p — 1 being the number of 
times the orbit visits U+ (recall that there is an extra 
symbol '1' coming from the gap). 

Suppose now that we introduce a small perturbation 
of the parameters and the initial condition: (xq, s, a) — > 
(x' ,e',a'). The perturbed orbit x' t — $*, a , (xq) is ini- 
tially at a distance Axo — x' — xq from the unperturbed 
one. Using the continuity of x' t on the parameters e, a 
and the initial condition, the distance Axt between the 
two orbits at time t can be made as small as we want by 
making (xq,e, a) sufficiently close to (x' ,e',a'). Because 
r is an open interval, if x q belongs to T, then so does x' q 
for a sufficiently small perturbation. Therefore we have 
established the crucial result: 

Theorem 3. If the gap orbit of $ eiQ is finite, then it is 
stable under a sufficiently small perturbation of parame- 
ters and initial condition. 

Thus there is a region in the parameter space (e, a), 
where the given rotation number is constant or mode- 
locked. An example of mode-locking e-regions for a fixed 
value of a is given in figure 1. The mode- locked region, 
corresponding to a given rotation number v — p/q, can be 
computed by noting that in order that an orbit through 
any point in T returns after q — 1 iterations we must have 
that $|~ 1 (r) C r. This condition may be rewritten as 
<I>f~ 2 (l — 2e) G r giving the inequalities 

7+ < <5>l: a 2 {l-2e) < 7 -, (31) 

since 

Vx er =s> xi = $ Eta (x ) = l-2e. 

We shall mention that the end points of ( |3l| ) are included 
because an orbit arriving at 7± after q — 2 iterations 
will also reach 1 — 2e after two more iterations since 
3 e ,a(/±(7±)) = *e,«(±a) = 1 - 2e (cf. ©) - these 
two orbits correspond to the kneading sequences of the 
local map jl(j|n|]. The closed e-interval region for which 
the velocity is mode-locked to a given rational is called a 
plateau. 

Equation (|3l|) ensures that the orbit starting at x\ = 
1 — 2e falls into the gap after q — 2 iterations and there- 
fore repeats itself after q — 1 iterations. The value of 
<I>f~ 2 (l — 2e) is determined by the coding sequence of the 
associated rotation number. So if we want to compute 
the mode-locked region for a given velocity v — p/q we 
have to determine the (e, a) region where ( pl|) is satisfied 
in correspondence to the code S(p/q) given by ( |30| ). For 
example, to compute the mode-locked region for the ve- 
locity v = 2/5 we first compute S(2/5) = (0,1,0,1,0), so 
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we know that the orbit visits the regions U-, U+, U- and 
r. Solving for $3 o (! _ 2s ) = f_ (/+ (/_(1 - &))) gives 

7+ < /-(/+(/-(l-2 £ ))) < 7 _. 

In figure 5 we show the mode-locking regions in the (e, a)- 
plane, computed via (j3l|), for some rational velocities. 
These regions, represented as shaded areas, are the so- 
called Arnold's tongues |l^Jlq| . It can be observed that 
as we approach the zero-gap line (dashed line) — given 
by 7 = in (|l9|) — the tongue size decreases, as the 
latter is related to the gap size. 



B. Lattice representation 

We now introduce a geometric representation of the 
rational velocities of the wavefront. This representation 
gives a straightforward method for computing the associ- 
ated coding sequence as well as some properties related to 
the mediant of two given consecutive rational velocities 
of a Farey sequence. 

We represent the rational velocity v = m/n as the 
point P(m/n) = (n — m,m) on the integral lattice Z 2 . 
Consider the ray OP emanating from the origin and pass- 
ing through P(m/n). It is obvious that distinct irre- 
ducible rational velocities correspond to distinct points 
and distinct rays. The ray OP will intersect the lattice 
grid infinitely many times. We drop the first intersec- 
tion in O and take all the remaining ones. Next, write 
a '0' every time the ray OP crosses a vertical grid line, 
write a '1' when it crosses a horizontal grid line and write 
'10' when its crosses through a lattice point (see example 
in figure 6). The resulting semi-infinite binary sequence 
can be shown to be equal to the symbolic coding of the 
velocity v = m/n given by (|3C|). 

We sketch the proof. Recall that a '0' in S(m/n) means 
one iteration of the CML without advance of the interface 
while a '1' means that the interface has advanced one 
site. Let us now consider the consequence of crossing 
vertically and horizontally the grid lines of Z 2 . Crossing 
a vertical line of the grid means transition from the point 
(n — m, m) to the point (n + 1 — m, to), i.e. from P(m/n) 
to P(mj (n+1)) — this corresponds to iterating the CML 
once without any advance of the interface (symbol '0'). 
Crossing a horizontal line means moving from the point 
(n — m, to) to the point (n — to, to + 1) = ((n+1) — (to + 
1), to + 1), i.e. from P(m/n) to P((m + l)/(n+l)), which 
corresponds to a shift of one site of the interface after 
one iteration (symbol '1'). And finally, crossing through 
a lattice point is equivalent to a shift of one place of the 
interface in two iterations (symbol '10'). 

~~ jHU of i 
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>rder r, %, is 



We recall that the Farey series 
defined to be the set of irreducible fractions, in ascending 
order, belonging to [0, 1] whose denominators are smaller 
than or equal to r. Suppose now that the velocities 



v = m/n and v' = m'/n' are two consecutive rationals of 
the Farey series 3m a x(n.n')- Then, using the lattice repre- 
sentation, it is not difficult to see that if m/n < m'/n' the 
coding sequence of the velocity v" = (to + m')/(n + n') 
corresponding to their mediant is given by 



S 



= s 



s 



(32) 



where the right hand side denotes the concatenation of 
two periodic sequences, defined as follows. If S a = 
(of] ■ ■ ■ , a P ) and Sb — (pi, ■ ■ ■ , b q ) are periodic sequences 
of period p and q, respectively, their concatenation S a Sb 
is the periodic sequence of period p + q defined as 



S a Sb — (ax, . . . ,a p ,bi, . . . , b q ). 



(33) 



The result (^2|) derives from the fact that the parallelo- 
gram (OPP"P'), with P = P(m/n),P' = P'(m'/n') and 
P" = P"(m"/n") = P"({m + m')/(n + n')), does not 
contain any lattice point since m/n and m'/n' are con- 
secutive numbers of a Farey series [ p0[ . 

An example of this construction is shown in figure 6 
where the velocities v = 1/3 and v' = 1/2 give the medi- 
ant v" = 2/5. The sequence for v = 1/3 is (010) (vertical 
cross + lattice point) and the sequence for v' = 1/2 is 
(10) (just a lattice point) and the resulting sequence for 
the mediant (01010) is obtained by the concatenation 
(010) (TO). 



C. Unimodular transformations and envelopes 

The structure of the parameter space of the map- 
ping F eja at the boundary of a tongue can be described 
analytically by means of envelopes. These are func- 
tions characterizing the structure of sequences of adja- 
cent tongues. The study of envelopes proceeds in two 
stages. Firstly we derive upper and lower envelopes 
for the so-called first order plateaus (i.e. v — 1/n and 
v = (n — l)/n, n — 1,2,.. .). Then with the use of uni- 
modular transformations we find envelopes for any family 
of plateaus. 

The zeroth order plateaus are defined as the two 
plateaus given by ( p4| ) with velocities {0/1, 1/1}. In or- 
der to construct any order of plateaus, the k-th order 
plateaus, we take any two consecutive elements of the pre- 
vious order k— 1 denoted by v a — m a /n a and Vb — mb/nb 
(v a < Vb). Then the next order of plateaus is defined as 
the following infinite increasing sequence of mediants 

3m a + TOh 2m a + m a + to^ m a + 2m,;, 
3n a + n b ' 2n Q + n b n a + n b n a + 2n b 



(34) 



The first order plateaus are thus given by 

1112 3 

''4'3'2'3'4'' 
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They form an unique family. On the other hand, between 
any two consecutive plateaus of a given order there is a 
new family of plateaus of the next order and so on. This 
construction is similar to that of the Farey series Jl9| . 
In particular, any two consecutive plateaus p/q and j//g' 
(p/q < p'/q') of the same order and family are consecu- 
tive fractions of a Farey series (i.e. qp' — pq' = 1). 

From the definition of the sequence of plateaus ( ]34| ) 
and by applying repetitively the sequence concatenation 
( |32| ) it is possible to find the sequences associated to a 
given family of plateaus: 

{ . . . , s a s b , s a s b , S a s b , s a , S a !?(,,... y , 

where we used the shorthand notation S a = S(m a /n a ) 
and Sb = S(mb/rib)- It will be useful to distinguish the 
left S a and right Sb subfamilies of sequences defined by 



S Q — {■ ■ ■ , S^S b , S^Sbi S a Sb} 
S& = {S a S b , S a Sfo, S a Sfo 7 . . .} . 



(35) 



Note that the sequence S a S b , corresponding to the medi- 
ant of v a and Vb, appears in both families. The elements 
of the subfamilies S Q (Sb) tend, to the left (right), to S a 
(Sb). Thus the subfamilies S (Sb) correspond to veloc- 
ities going from the mediant of v a and Vb to the plateau 

V a (V b ). 

From now on we consider only S a . All the results 
generalize to S& just by interchanging the subscripts a 
and b. The elements of S n are of the form 



SeS a => S = S2S b n=l,2,. 

which correspond to the velocities 
m a + m b 



v(S:S b ) 



rib 



1,2, 



(36) 



(37) 



We recall that the transformations (|37j) are unimodular 
since we always deal with consecutive fractions of Farey 
series. 

The sequences S a and Sb correspond to a particular 
combination of /_ and / + . Since f± are linear, any such 
combination will again give a linear function. We de- 
note by fs a (fs b ) the linear function resulting from the 
combination of f± specified by the sequence S a (S b ): 



fs a (x) = a a x + (3 a 
fs b (x) =a b x + f3 b . 

Composing fg a with itself n times gives 

/*>) a +c 



(38) 



(39) 



where £* = /3 a /(l — a a ) is the fixed point of fg a . Com- 
bining ( |38| ) and ( |39| ) in (36) we obtain the iterates of 
xo = 1 — 2e under & E:a for the coding sequences S%Sb 



Replacing @ in ^ we obtain the mode-locking region 
for a given series of plateaus 

/io(7+) < O6(<£(zo -£)+£)+ & < /io(7-), 

(41) 

where we have taken into account the fact that the two 
final symbols '10' of any sequence in S a come from the 
gap r and therefore they have to be compensated by 
replacing j± in @) by /i(/o(7±)) = /io(7±) ■ 

As we mentioned earlier, the plateaus include their end 
points. These are obtained from solving the inequalities 
( p| ) for the extremes /io(7±)- Therefore, all the end 
points of a given family of plateaus are given by 



a b (aZ(x o -£) + £) + 0b = Ao(l±), 



(42) 



where the positive (negative) sign gives the left (right) 
end points of the plateaus. As we vary n — 0, 1, . . . in 
( [i"2| ) we browse all the extremes of the series of plateaus 
going from the mediant of v a and Vb to v a . 
Solving (|42|) for n yields 



n±(e,a) 



1 



In < 



In 



/lo(7±) ~ Pb~ UbCa 



Finally, replacing the value of n in (| 
the continuous functions 



by n±(s,a) gives 



v±(e,a) 



n±(e, a) m a + mt 
n±(e,a)n a + n b 



(43) 



fs S sM = f Sb (fl(xo)). 



(40) 



that passes through all the left (v + ) and right («_) 
boundaries of the plateaus. Therefore, v± gives the up- 
per (+) and lower (— ) envelopes of the series of plateaus 

In figure 7a we display the first order plateaus for a — 
0.4. We only plotted the left family (v < 1/2) since the 
right family (v > 1/2) is symmetrical from (0). An 
example of a second order family, for the same value of 
a, is shown in figure 7b where we display the lower and 
upper envelopes for the left and right families of second 
order plateaus between v = 1/2 and v = 1/3. 

With the method described above, it is possible to find 
the envelopes of any sequence of high order plateaus via 
unimodular transformations. The self-similarity struc- 
ture of the velocity function, easily observable in figures 
7a and 7b, is controlled by modular transformations d37j ) 
and envelopes. Since it is possible to find a periodic 
symbolic sequence for any given rational velocity arising 
from an eventually periodic orbit of $ e a through 1 — 2e 
then it is possible to find a non-empty e interval where 
the velocity is mode- locked (Theorem 3). Therefore the 
graph of the velocity v(s,a) as a function of e is in fact 
a Devil's staircase p^j2l| ], i.e. a fractal staircase. This 
mode-locking phenomenon is very similar to the one ob- 
served for a uniform rotation in a subcritical perturbed 
circle map p4l. 
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If, on the other hand, the symbolic sequence is chosen 
to be non periodic, the associated velocity v is irrational. 
The continuity of the function v(e) guarantees that there 
exists at least one value of e such that v(e) = v. Wc 
believe that it is possible to prove that there is only one 
such e. The idea is to take a sequence of rational approxi- 
mants of v; the corresponding upper and lower envelopes 
should approach each other, as the approximants tend to 
is, and coincide in the limit to give a unique e such that 
v(e) = v. 



D. Convergence to minimal states 

In this chapter we have described the structure of the 
parameter space associated with minimal states. Here we 
generalize the notion of minimal state and briefly com- 
plete our description of the parameter space. 

A minimal state — a step state of size or 1 — is stable 
under iteration when the gap is non-negative. When the 
gap is negative, the stable configuration of a step state 
is found to have size assuming only two values, N and 
N — 1. Such states will be called minimal A-states. 

Numerical experimentation reveals that the integer iV 
does not depend on the initial condition. Any initial 
configuration we have tried in the interface (increasing, 
decreasing, oscillatory, random with small and large size, 
etc . . . ) converged, after a transient (depending on the 
initial conditions) to a minimal A-state, with the value of 
N depending only on e and a. Accordingly, the param- 
eter space is partitioned into layers where the value of 
N is constant, as depicted in figure 8. The boundary be- 
tween the first layer (minimal 1-state) and the second one 
(minimal 2-state) is the zero-gap curve, given by 7 = 
in ( |l9|) . Above that curve the gap is negative, and there 
appears to be an infinite family of layers labeled by the 
value of N, separated by boundary curves which may not 
be differentiable. 

The dynamical analysis for such minimal A-states may 
be carried out by reducing the dynamics of the whole 
lattice to a A-dimensional auxiliary map |2^| analogous 
to the one obtained previously for minimal 1-states. 



VI. GENERALIZATIONS AND CONCLUSIONS 

We have described a mechanism that allows travelling 
interfaces in a one-way CML (J5J) for a special family of 
piece- wise linear maps (^). We now give some ideas on 
the generalization of these results to a broader range of 
functions and couplings. 

First of all, the superattractiveness of the points x*_ 
and x* + is not necessary. It suffices that the map / be 



a contraction mapping at x±. This is the case, for in- 
stance, if f"(x) < for x*_ < x < x* and f"(x) > 
for x* < x < x* + , where x* is the unstable fixed point in 
between x*_ and x* + - However, without superattractive- 
ness, the travelling wave will be less stable to perturba- 
tions. Then, if we are dealing with a medium involving 
random fluctuations we would like the map to be suffi- 
ciently contracting for a travelling interface to survive. 
On the other hand, the lack of superstability of x± intro- 
duces more sites in the interface since the convergence to 
x*± is slower. Moreover, in this case the number of sites 
contained in the open interval (x*_ ,x* + ) is not necessar- 
ily finite (take the example of an exponentially localized 
state). For this reason the reduction of the dynamics in 
this case does not lead to any simplification. Neverthe- 
less, if the shape of the front remains unchanged during 
the evolution, it is possible to concentrate our attention 
to a single site in the interface since the remaining sites 
are determined by the shape of the front. We could, for 
example, consider the site that is closest to the unstable 
point x* and from it reconstruct a one-dimensional aux- 
iliary map that accounts for the whole dynamics, in the 
same way that <£> eia accounts for the whole dynamics of 
minimal states. A more detailed study will be given in 
II- 

Secondly, the discontinuity of the derivative of the local 
map is not necessary for mode-locking. We have observed 
mode-locking for analytic maps as well. In figure 9 a win- 
dow of the velocity as a function of e is shown for some 
polynomials maps (of odd degrees from 5 to 13). The 
plateaus v = 1/2 and v = 1/3 are clearly visible. The 
origin of mode-locking is not related to non-analyticity. 
It ultimately lies with the asymptotic (t — > 00) discon- 
tinuity from the initial conditions associated with the 
repelling fixed point x* between the two stable points, 
which separates two dynamical behaviours. 

In this respect the value of the derivative at the repeller 
can be interpreted as a measure of this discontinuity. It 
is therefore expected that, for fixed value of the deriva- 
tive at the attractors, the mode-locked plateaus become 
more pronounced as the derivative at the repeller point 
increases. In figure 9 the plateaus v = 1/2 and v = 1/3 
are almost imperceptible for the polynomial of degree 5 
that has a derivative at the repelling point ~1.9, mean- 
while they become more appreciable for the polynomials 
of higher degree (7, 9, 11 and 13) where the derivative 
at the repelling point gets larger (~2.2, 2.5, 2.7 and 2.9 
respectively). 

Finally, the travelling wave property arises from one- 
way coupling in the CML and the direction of the wave is 
given by that of the coupling. By comparison, a front on 
a diffusive CML ([|) with a symmetric map at each site is 
going to remain stationary in the long term. For the front 
to advance one needs asymmetry, either on the coupling 
or in the local map (with respect to the repeller point x*) 
which generates a bias between competing attractors. 
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In order to illustrate this phenomenon let us take the 
logistic map f^(x) = fxx(l — x). The logistic map itself 
does not fulfill the requirements (two attractive points 
with a repeller in between) but consider its second iter- 



ate ft 



1 + V5 the map- 



At the parameter value p 
ping / 2 has two superattractive fixed points x*_ = 1/2 

and x* + = 1/2 + 3 — \/5/4) (corresponding to a 
superattractive 2-cycle for /) separated by a repeller 
x* = 1 — s/2 Vo — \/5/4. Consider iterating the diffu- 
sive CML with the map at each site. For an initial 
step state we obtain a travelling interface where the ve- 
locity again behaves as in the examples shown above — it 
has a transition from v(e) = to vie) > at e c ~ 0.198 
— , see figure 10. In this case the mode-locking is not 
apparent at first glance. However, on amplification, the 
plateaus are shown to occur. For example, the plateaus 
corresponding to v = 1/7 (right) and v = 1/8 (left) are 
shown in the enlargements in figure 10. 

In conclusion, the mode-locking of the velocity of 
the travelling interface, appears to be an universal phe- 
nomenon in coupled map lattices. It implies that waves 
travelling with rational velocities are stable under both 
spatial and parametric perturbations. 

Appendix A 

In this appendix we provide various statements and 
proofs concerning the image of localized states. 

Lemma A.l. Let Xt be a localized state of a one-way 
CML, let x± = \im l ^ ±OD x t (i), and let M(X t+ i) ^ 0. 
Then, if f is bounded and continuous at x±, X t +i is a 
localized state. If, in addition, X t is exponentially local- 
ized, so is X t +\. 

Proof. We begin considering the case i — * oo. Because X t 
is localized, x+ is finite. Moreover, since / is continuous 
at x+ , then there exists a constant p > and an integer N 
such that, for alii > N we have \f(xt(i+l)) — f(xt(i))\ < 
p\xt(i + 1) — Xt(i)\, whence 

Ax t+1 (i) = |(1 - e) f(x t (i + l))+e f(x t (i)) 
-(l-e)f(x t (i))-ef(x t (i-l))\ 

< (l-e)\f(x t (i + l))-f(x t (i))\ 
+ e\f(x t (i))-f(x t {i-l))\ 

< p(l - e) \x t (i + 1) - x t (i)\ 
+ pe \x t (i) - x t (i - 1)| 

= p(l-e)Ax t (i) + peAx t (i-l), i > N. 

(A.l) 

Then, from absolute convergence 

Y,Zn+i Ax t+i(i) < p(l-e) ££U+i Ax t (i) 
+ P £ E~ jv+i Ax *(* _ !) < °°. 



A similar inequality holds for i — > — oo. From the above 
and the fact that any finite sum of Axt+i(i) is bounded 
(because / is bounded) it follows that M t +i is finite (and 
non-zero, by hypothesis), so that the probability pt+i 
exists. Finally, multiplying both sides of the inequality 
(A.l) by i and by i 2 , respectively, and summing over 



i > N shows that (it+i and cr t + 1 are finite, whence X t +i 
is localized. 

If Xt is exponentially localized, then there exist posi- 
tive constants c and n < 1 for which 



Ax t (i) < ck\ 



(A.2) 



Combining the above with (A.l) we obtain 

Ax t+1 (i) < K l cp ((1 - e) + ^) , i>N. (A.3) 

A similar estimate holds for negative i, showing that 
Xt+\ is exponentially localized (possibly with a larger 
constant c). □ 

From Lemma A.l and the fact that the mass of a step 
state is bounded from (cf. (||)) we obtain immediately 

Corollary A.2. Let X t be a step state of a one-way 
CML, and let f be continuous at x* ± and bounded. Then 
X t +\ is a step state. If, in addition, X t is exponentially 
localized, so is X t+ i . 

We have then establish the invariance of the localized 
and exponentially localized states under the dynamics 
of the CML and we are ready now to give a proof for 
Theorem 1. 

Theorem 1. Let Xo be an exponentially localized step 
state of a one-way CML. Let the local map f be bounded, 
and let f be a contraction mapping in a neighbourhood 
of the fixed points ([?]). Then, for all sufficiently small e, 
v(e) = and v (1 — e) = 1, independently of Xo. 

Proof. We deal with the case of small coupling first, 
and positive i. From Corollary A.2 we know that X t is 
an exponentially localized step state for all t > 0. By 
assumption / is a contraction mapping in some domain 
\x — x* + 1 < r. Let N' be such that for all i > N' we have 
\x t {i) -x* + \ < r/2, so that Ax t (i) < r. Then, for i > N', 
the constant p appearing in ( A.l ) can be chosen to be 
smaller than 1. 

Let TV" be such that Ax t {i) < 1 for all i > N" , and 
let N = max(N',N"). Then since X t is ex pone ntially 
localized, we can choose n so that the bo und (A.2 ) holds 
for i > N with c = 1. Letting c = 1 in ( A.3 ) it follows 
that 



Ax t+l {i) < k 1 
provided that 



i > N 



(A.4) 



K 1 
£ < - 



P 1 
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which is satisfied for sufficiently small e, since the right 
hand side is positive. 

It remains to consider the case i — N. Because / is 
bounded, the quantity 

A/ = sup x f(x) - M x f(x) 



is finite. We have, in place of ( A.l) 

Ax t+1 {N)<p{l-e)Ax t (N)+eAf 
giving, for all i > 

Ax t+1 (N) < a t+1 n N + eAf 1 



v*+l 



1 — a 

where a — p{\ — e) < 1. The right-hand side of the 
above inequality can be made smaller than k n for all t, 
prov ided that e is sufficiently small. Thus the inequality 
( \KJ${) holds also for i = N. 

Let i be negative. Assuming that Axt{i) < kT % for all 
sufficiently large (negative) i, and proceeding as above, 
we find that the bound Axt(i) < kT % for i < —N implies 
that 



Ax t+ i(N) < k l p(l-e + SK,)<n l ; 



i < -N 



for all £ (all constants have been redefined). The case 
i = —N gives the recursive inequality 



Ax t+1 (-N) < (1 - e) A/ + pe/s- 



-N+l 



which, as before, yields Axt+i(— N) < k~ n for suffi- 
ciently small e. 

The above induction establishes the time-independent 
bound 



Ax t {i) <n ll{ , 



\i\>N, t>0. 



for a suitable choice of k and N, and for all sufficiently 
small e. In this e-range, from the boundedness of / we 
conclude that \pt\ and at are bounded from above for all 
times. Thus v(e) — 0, and X t remains a step state also 
in the limit t — > oo. □ 

Appendix B 

Proof of Theorem 2. 

(i) The dynamics of the i-th site is given by Xt+k(i) — 
f-(xt(i)) since it is coupled to the (i — l)-th site whose 
value is — 1. For the dynamics of the (i + l)-th site, 
since e < e c and a < f+(x) for any x € U, we have 
that Xt+k{i + 1) € S+. Therefore at time t + k the state 
of the lattice will be equivalent to X t +k = [f-{xt),i] 
which depends exclusively on the initial condition Xt(i). 
If e < e c the derivative of the map f'_(x) = (1 — e)/a > 1 
and then /_ has a repeller at Therefore f_(xt(i)) 
will decrease (increase), for Xt(i) < X- (xt(i) > xJ) until 
it reaches S- (S+) when the state of the lattice is P(i) 



(P(i — 1)). For the marginal case x t (i) — x_ the state of 
the lattice X t +k — [f-{ x -),i] = [f-(x-),i] is fixed and 
unstable since X- is an unstable fixed point of /_. 

(ii) Using a similar reasoning, if e > 1 — e c , x + is an 
unstable fixed point of / + and the state of the lattice at 
time t + k is, X t +k = [fX{ x t),i + k]. Therefore the lattice 
will reach the state P(i + k) (P(i + k— 1)) for x t (i) < x+ 
(xt(i) > x + ). In the marginal case x t (i) = x+, the 
state of the lattice is X t+ k — G k i + kf) — 

G k ([f+(x+], i]), which is fixed and unstable under GoF. 
□ 
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Figure 1: Travelling velocity v(e) of the interface, for a 
piece-wise linear one-way couple map lattice as a function 
of the coupling parameter e, computed numerically for dif- 
ferent values of the local map parameter a. The velocity has 
a staircase-like behaviour symmetric with respect to the point 
(e,v(e)) = (1/2,1/2). 
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Figure 2: Piece-wise linear local map f a for the one-way 
CML. 
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Figure 3: The auxiliary map $ e a which provides the dynam- 
ics of the interface. 
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Figure 4: Reduction of the degrees of freedom for minimal 
states of a one-way CML. The originally infinite-dimensional 
CML is reduced to a 1-dimensional system. The only degree 
of freedom is the value of the interface site, its dynamics is 
being given by the auxiliary map 




Figure 5: Arnold's tongues in the (e, a)-plane corresponding 
to rational velocities of the travelling interface for the one-way 
coupled map lattice F e a . 




Figure 6: Integral lattice representation of the rational ve- 
locities. Every rational velocity m/n corresponds to the point 
P(m/n) = (n — m, m) on the lattice. Its symbolic sequence is 
then constructed by the crossings of OP with the grid lines: a 
'0' corresponds to crossing a vertical line, a '1' a horizontal line 
and a '10' corresponds to crossing a lattice point. Concatenat- 
ing sequences for two different velocities m/n and m' /n' cor- 
responds to the vector addition of OP (m/n) and OP(m' / n 1 ), 
and the resulting velocity is the mediant (m + m')/(n + n'). 
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Figure 7: a) Upper and lower envelopes of the first order 
plateaus for a = 0.4. b) Upper and lower envelopes of the 
second order plateaus between v(e) = 1/2 and v{e) = 1/3. 
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Figure 8: Convergence layers for minimal states in the piece- 
wise liner CML. A TV-minimal state is a state with N or N — 
1 sites on the interface. The curve separating the layers 1- 
minimal state and 2-minimal state is called the zero-gap curve 

( 7 M) = o). 
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Figure 9: Mode-locking of the velocity for a CML 
with smooth polynomial local maps P n . The polynomials 
P5, P7, P9, Pn and Pi 3 have degrees 5, 7, 9, 11 and 13 respec- 
tively. The plateaus v = 1/3 and v = 1/4 become increasingly 
visible with the degree of the polynomials corresponding to 
increasing the derivative at the origin. 




Figure 10: Mode-locking of the velocity of the travelling in- 
terface in a diffusive CML. The local map is the second iterate 
of the logistic map f^(x) = fix{l — x) with \i = 1 + and 
the interface is supported by its two superstable points. The 
velocity seems to have a smooth graph, but the mode-locking 
plateaus corresponding to v = 1/7 and v = 1/8 become evi- 
dent under magnification. 



